Phases of a conserved mass model of aggregation with fragmentation at fixed sites 
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To study the effect of quenched disorder in a class of reaction-diffusion systems, we introduce a 
conserved mass model of diffusion and aggregation in which fragmentation occurs only at certain 
fixed sites. On most sites, the mass moves as a whole to a nearest neighbour while it leaves the 
fixed sites only as a single monomer (i.e. chips off). Once the mass leaves any site, it coalesces with 
the mass present on its neighbour. The total mass of the system is conserved in this model as is 
evident from the rules. We study in detail the effect of a single chipping site on the steady state in 
arbitrary dimensions, with and without bias. In the thermodynamic limit, the system can exist in 
one of the following three phases - (a) Pinned Aggregate (PA) phase in which an infinite aggregate 
(i.e. an aggregate with mass proportional to the volume of the system) appears at the chipping 
site with probability one but not in the bulk, (b) Unpinned Aggregate (UA) phase in which both 
the chipping site and the bulk can support an infinite aggregate simultaneously, (c) Non Aggregate 
O ' (NA) phase in which the system cannot sustain a cluster of mass proportional to the volume at 

all. The steady state of the system depends on the dimension and drive. A sitewise inhomogeneous 
Mean Field Theory predicts that the system exists in the UA phase in all cases. We tested this 
prediction by Monte Carlo simulations in Id and 2d and found that it is true except in Id, biased 
case. In the latter case, there is a phase transition from the NA phase to the PA phase as the mass 
■ density is increased. We identify the critical point exactly and calculate the mass distribution in 

the PA phase. The NA phase and the critical point are studied by Monte-Carlo simulations and 
using scaling arguments. A variant of the above aggregation model is also considered in which total 
particle number is conserved and chipping occurs at a fixed site, but the particles do not interact 
I . with each other at other sites. This model is solved exactly by mapping it to a Zero Range Process. 

' O ■ With increasing density, it exhibits a phase transition from the NA phase to the PA phase in all 

dimensions, irrespective of bias. The model is also solved with an extensive number of chipping sites 
with random chipping rates and we argue that the solution describes qualitatively the behaviour of 
the aggregation model with extensive disorder. 
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I. INTRODUCTION 



Reaction-Diffusion systems form an important class of nonequilibrium systems whose dynamics and steady state 
depend on various factors such as the nature of the reaction (aggregation, annihilation, birth, fragmentation), number 
or type of reactants involved (single or multi-species), velocity of reactants (ballistic or diffusion controlled) and 
presence of external input (injection) Q-j^] . An interesting class with wide ranging applications involves the elementary 
moves of aggregation (coalescence on contact) and fragmentation (break-up of clusters of masses), besides diffusion. A 
I ■ number of analytical results including the occurrence of nonequilibrium phase transitions have been obtained for such 
systems with translationally invariant geometries A natural question arises - what is the effect of quenched 

disorder on the possible phases of such systems and the transitions between them? We may anticipate interesting 
, effects, as quenched disorder is known to strongly influence the character of the steady state in other nonequilibrium 
systems P-|t1| . For instance, the steady state of driven diffusive systems on a one-dimensional lattice with bias shows 
phase separation in the presence of even a single defect [jlO 15-1^1. 



To see what effect quenched disorder might have in a system with diffusion and aggregation, consider the process of 
polymerisation in a random medium with traps at certain fixed sites in which the polymer can get stuck. Aggregation 
occurs when two diffusing chains meet and coalesce; the reduced mobility of the aggregates at the trap sites promotes 
the formation of large, localised aggregates at such sites. If these traps are not perfect and allow monomers to detach 
and leave, there is also a possibility of formation of large mobile aggregates in the bulk. 

To elucidate under what circumstances which types of aggregates may form as a result of these physical effects, 
namely, diffusion-aggregation, chipping and trapping, we consider a simple, reduced model in which the mass leaves 
as a single aggregate to a nearest neighbour from all sites except at certain, fixed sites from which it is allowed to leave 
only by single monomer dissociation. We refer to the loss of unit mass as chipping and these special sites as chippers. 
The total mass in the system is conserved as is evident from the dynamical rules described above. The quenched 
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character of the chipper sites is important in determining the steady state of this model which is quite different from 
that in the translationally invariant uniform-chipping model where chipping can occur at every site jTj. 

As a first step towards understanding such a spatially inhomogeneous system, we study in detail the case when 
only a single chipper is present. Despite the simplicity of the model, we find that that the presence of a chipper gives 
rise to interesting steady states. We study the model both analytically and numerically in arbitrary dimensions, both 
with and without a global bias which sets up an overall mass current. The constraint of conservation of total mass 
in the model plays an important role in determining the steady state. We find that the system may exist in three 
possible phases. Depending on the dimension and the presence of drive, the system either stays in one of the phases 
or else make transitions from one phase to the another as the parameters are varied. These phases are characterised 
by the presence or absence of aggregates with mass proportional to the volume of the system. In the thermodynamic 
limit, the mass of such an aggregate diverges, so we refer to it as an infinite aggregate. The three phases are described 
below: 

Pinned Aggregate phase (PA phase): In this phase, with probability one, an infinite aggregate occurs only at 
the chipper site and not in the bulk. This infinite aggregate acts as a spatially localised particle bath for rest of the 
system. 

Unpinned Aggregate phase (UA phase): In this phase, a localised infinite aggregate can exist at the chipper 
site together with mobile infinite aggregates in the bulk. Interestingly, both types of infinite aggregate can exist 
simultaneously. In this sense, this state is different from the aggregate phases in translationally invariant systems 
such as the no-chipper limit of this model or the uniform chipping model [Q , in which only one mobile infinite aggregate 
occurs at a time. 

Non Aggregate phase (NA phase): This phase is characterised by the absence of an infinite aggregate anywhere 
in the system. The mass is spread out all over the system in clusters, each of which has a vanishing fraction of the 
total mass in the thermodynamic limit. 

We analysed the system within a Mean Field Theory (MFT) allowing for the spatial dependence in the mass 
distributions. We find that the MFT predicts the existence of only the UA phase. On comparing this result with 
Monte Carlo simulations in d = 1 and d = 2, we find qualitative agreement except in the lei, biased case. In this 
exceptional case, there is a phase transition from the NA phase to the PA phase as the density is increased. This 
case is different from the rest due to an interesting interplay of the ballistic scale of motion and the diffusive scale of 
coalescence in one dimension. 

In the model described above, an infinite aggregate is formed due to both interactions (coalescence of diffusing 
particles) and the presence of an inhomogeneity in the form of the chipper. It is useful to contrast the behaviour of 
this model with that of a model of noninteracting particles in which the aggregate is formed solely due to the presence 



of disorder. We solve this model in the presence of a single defect by mapping it to a Zero Range Process |L4 18 
and show that it exhibits a phase transition from the NA phase to the PA phase in all dimensions, irrespective of 
bias. This is in contrast to the aggregation model described earlier which does not show a phase transition in higher 
dimensions in the presence of single chipper. 

Although we mainly discuss the case of a single chipper in this paper, we also discuss what happens in the presence 
of extensive disorder in both the models. We argue that in the presence of an extensive number of chippers, the 
interaction effects (i.e. coalescence) can be ignored in the aggregation model on large enough length and time scales 
and it behaves like the free-particle model. The latter model can be solved in the presence of extensive disorder and 
shows a phase transition in all dimensions for all bias from the PA phase to the NA phase as the density is decreased. 

The remainder of the paper is organised as follows. We define the single-chipper aggregation model in Section || and 
discuss the possible phases on the basis of the conservation law. We analyse the system within a MFT in Section |l| 
and show that it predicts the occurrence of the UA phase. We also discuss the numerical results which support this 
broad conclusion in several cases. In Section [w| we turn to the exceptional case namely the one dimensional, biased 
case. We present analytical results in the PA phase and numerical results in the NA phase and at the critical point. 
In Section we present the solution of the single-chipper free particle model. We also discuss the likely behaviour 
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of the system with extensive disorder in Section Finally we summarise our results in Section 

II. SINGLE-CHIPPER AGGREGATION MODEL AND ITS PHASES 

A. The Model 

Our model is defined on a <i-dimensional hypercubic lattice of length L with periodic boundary conditions, and is 
studied both in the presence and absence of bias. The system evolves via following rules: at any site except the origin, 
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FIG. 1. Illustration of the aggregation and chipping moves in the totally biased case in one dimension. In an elementary 
move, only a monomer (shaded) can leave the chipper site (located at the origin) at a rate w. At sites other than the origin, 
the mass leaves as a whole (shaded) at a rate 1. Once the mass moves from any site to its neighbouring site, it coalesces with 
the mass already present on that site. 



all the mass at that site moves as a whole to a nearest neighbour at a rate 1, while at the origin only a single monomer 
chips off at a rate w leaving rest of the mass behind. Once the mass moves from any site to its neighbouring site, it 
simply coalesces with the mass present on that site (see Fig. |l|). Evidently, the total mass M = pL d of the system is 
conserved, where p is the mass density. The evolution of the system can be described by the time evolution equation 
for the probability Pfe(m, t) that there is mass m at site k at time t. As the moves at the origin are different from that 
at the rest of the system, the mass distribution at the origin (denoted by 0) and its nearest neighbours which receive 
mass from the origin (denoted by set A) obey equations different from those obeyed by rest of the system (set B). 
The set of time evolution equations obeyed by the sites in the bulk (set B) are given by 



dP k (m,t) 
dt 

dP k (0,t) 
dt 



y^[^ Pk,k+s(m- m',m',t) - ^ P kl k+s(m,m',t) - P k (m,t)] , k e B (la) 

8 m' — l m'^0 

^P^tl-^P^mVll.^fl , (lb) 

5 m/0 m'/0 



where P k k+ g(m, m', t) is the joint probability that site k and its nearest neighbour k + 5 have mass m and m' 
respectively. The index S runs over the 2d nearest neighbours in the unbiased case and d forward neighbours for 
totally biased case. The convolution term in Eq.([la|) is the gain term through which the deficit mass is supplied to 



site k via diffusion-aggregation moves from its nearest neighbours. The gain term for zero mass in Eq.(lb) results 
from hopping the mass from site k to one of its nearest neighbours. The system can get out of the configuration in 
which the site k has mass m (including zero) when (a) the mass hops into site k from its nearest neighbours and (b) 
the nonzero mass at site k hops to its nearest neighbours. 

For the sites in the immediate neighbourhood of the origin (set A), the time evolution equation for P k (m,t) is 
similar to that for set B except that the contribution of the origin needs to be taken into account separately. We have 



8P k (m,t) 



m 

ni = X!' X! p k,k+s(m- m\m!>t) +w 22 ^fe,o( m ~ ^^t) - }/ ^ P k . k+ s(m,m\t) 

8 m' — l m'^0 8 m'^0 

-w Pkfl{m,m',t) -Y^Pk(m,t) , k € A (2a) 

m'=£0 5 
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dP k (0,t) 

at 



= J2J2 P ^ m ^-J2' Z) p k,k+s{0,m',t)-w PkA°> m ',t) , keA , (2b) 



where the primed sum denotes the sum over the nearest neighbours, excluding the origin. The contribution of the 
origin is taken care of by the terms with a coefficient w which accounts for the gain or loss in mass via chipping from 
the origin. 

At the origin itself, the set of evolution equations obeyed are given by 



dP (m,t) 
dt 

dP (0,t) 



<■)/ = /J?, Pofl+sjm - m', to', t) + wP (to + 1, t) - } j P ,o + g(m, to', t) - wP (m, *)] (3a) 

S in' — 1 rn'^O 

m = £>i , o(M)- Y, Po,o+s(0,m',t)] . (3b) 

S m'=£0 

In the gain term for nonzero mass, besides the convolution term, there is an extra term due to the possibility of 
chipping off one extra monomer with a rate w. The loss terms are similar to those discussed for the sites in set B. 

B. The Phases 

One can obtain useful information about the nature of the steady state from the following simple analysis. In the 
steady state, the mass current into and out of any site must be equal, i.e. 



(m fc ) = Y ( m k+s) , keB (4a) 
s s 

(m fc ) =ws + Y^ '( m k+s) , keA (4b) 
s s 

Y ws o = Y (mo+s) , (4c) 
s s 

where (m&) is the average mass at site k ^ and sq = 1 — Pq(0) is the probability that the origin is occupied. The 
solution of the above equations gives (OTfe) = wsq for all k ^ 0. Thus the average mass on a site in the bulk is uniform 
in spite of broken translational invariance. Since the total mass of the system is conserved, the average mass at the 
origin (Mo) must satisfy 

(M ) = P L d - ws (L d - I) . (5) 

Using this constraint equation, one can deduce some characteristics of the possible phases in the system as follows, 
bmcc the LHS in Eq.(§) is non-negative and sq is bounded above by one, the allowed domain for sq is constrained as 
shown in Fig. |[ Depending on the value of so (and hence (Mo)), the system can exist in one of the three distinct 
phases described below. 

PA phase: The system exists in this phase when sq is pinned to its maximum value 1 for all p > w (shown by the 
bold line in Fig. ||). From the constraint equation Eq.(||), one deduces that the average mass at the origin grows as 
L d giving rise to an infinite aggregate in the thermodynamic limit. Thus the Pinned Aggregate phase is defined as 
that in which an infinite aggregate is present at the origin with a probability one. 

UA phase: The system exists in this phase for all values of sq lying in the shaded region in Fig. ||. This phase 
is characterised by sq strictly less than 1, and not equal to p/w. As in the PA phase, the average mass at the origin 
grows as L d but Sq is not fixed and varies with p. Since sq < 1, the infinite aggregate is present at the origin only for 
a finite fraction of the time. For this reason, we refer to this phase as the Unpinned Aggregate phase. 

NA phase: This phase is characterised by sq — > p/w as L — * oo (shown by the solid line in Fig. ||). Using Eq.(||), 
one immediately obtains (Mo) /L d — > in the thermodynamic limit. Due to the absence of an aggregate with mass of 
0(L d ) at the origin, we call this the Non Aggregate phase. 

A more detailed study of the system shows that depending on the dimension and drive, the system can either stay 
in one of the three phases or else make transitions from one phase to another as p is varied, keeping w fixed. We have 
solved for (Mo) and sq in d dimensions within a sitewise inhomogeneous Mean Field Theory (Fig. ||) which predicts 
that the system exists only in the UA phase in all dimensions, regardless of the bias. This prediction was tested 
numerically in several cases and seen to be qualitatively correct except in the Id, biased case. 
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FIG. 2. Plot of so vs. p to show the three possible phases in the steady state of the aggregation model in the thermodynamic 
limit. The system exists in (a) Pinned Aggregate phase when so is one (bold horizontal line) (b) Unpinned Aggregate phase 
when so < 1 and not equal to p/w (grey striped region) (c) Non Aggregate phase when so is equal to p/w (dashed line). The 
smooth curve shows the MFT solution for so = p/(p + w) in the aggregation model. 



In the latter case, there is a phase transition from the NA phase to the PA phase as one crosses the critical line p = w 
either by increasing p with w held fixed or by decreasing w for fixed p. 



III. UNPINNED AGGREGATE PHASE: MEAN FIELD THEORY AND MONTE-CARLO SIMULATIONS 

In this Section, we analyse the aggregation model within the mean field approximation. We implement this ap- 
proximation by replacing the joint probability distribution P/^/j+^m, m' ,t) by the product Pk(m, t) Pk+s (m' , t) of the 
single site probability distributions. We begin with the MFT analysis in arbitrary dimensions in Section [III A|. We 



will see that MFT predicts only the existence of the UA phase in all dimensions, irrespective of bias. In Section IIIB, 
we study the UA phase in one dimension in more detail. We discuss some limitations of this MFT and also present 
some numerical results. We close this Section with a qualitative discussion of the UA phase. 



A. MFT in arbitrary dimensions 

In this Section, we solve for (Mo) and sq in arbitrary dimensions and for all bias within MFT and show that it 
predicts the occurrence of only the UA phase in all cases. 

In the steady state, the mass distribution at any site k is independent of time and is determined by setting the LHS 
of Eqs.(0)-@ equal to zero. Further we replace Pfc,k+{(ni, m') by Pfc(m) Pk+s(m') in Eqs.(f[a|)-(|3b|) to obtain the 



mean field equations for Pk(m). To study these equations, we define Qk{z) = J2 m =o Pk(jn)z m , the Laplace transform 
of Pfc(m) with respect to mass, and find that it obeys the following set of equations, 

^Q fe = ^[l + Q fe (Q fc+5 -l)] , keB (6a) 
s s 

Qk = vTT n 777] TT ' k e A ( 6b ) 

1 + ws {l -z)-\, s '{Qk+s - 1) 

y ° £>(i-*) + *(Qo+*-i)] ■ 1 ' 
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Since Pk(m) is normalised to one, Qfc(l) = 1 for all k in the above equations. The average mass in the bulk 
(mjfe) = Q' k (z)\ z= i still obeys Eqs.(f4a|)-(f4c|) which gives uniform mass in the bulk (m) = wsq, and Eq.(^) still holds. 

Although we could not solve the above set of nonlinear coupled equations, we were able to compute two quantities 
of primary interest, namely, (Mo) and sq. We begin by observing that the average mass at the origin can be written 
in terms of the mean-squared mass at its nearest neighbours, 

/W i T,s ( 2ws q + ( m l+s) - H) ™ 

EaM 1 - so) 

where we have used that (M ) = Q' (z)\ z =i and (mfc) — Q'l(z)\ z =\ + (m), k ^ 0. 

The mean-squared mass at sites other than the origin obeys the following set of linear equations, 

E< m fc>=E(Vfc+*>+ 2 <™> 2 ) > k ^ B ( 8a ) 

s s 

E ( m D = E ' (< m *+«> + 2 (™) 2 ) + 2 <™ 2 > + (m) , fc e ^ . (8b) 

8 8 

It suffices to calculate J^s ( TO o+<5) m or der to obtain (Mo). Adding the above L d — 1 equations , one obtains 

E ( m w) = E ( 2 <™> V - !) + (»»>) . ( 9 ) 

s s 

which further yields 

{Mo) = ™o + (rn)\L«-l) =pLd _ d _ 
w(l - s ) 

where the mass conservation equation Eq.([5]) has been used in the last identity. Solving for so and (Mo) in terms of 
p, w and L, we find 

x 1 (1 - w) + (p + w)L d v ' 

pL d 

s ° = 71 \ I / I — ■ ( llb ) 

(1 — w) + (p + w)L a 
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FIG. 4. Data collapse for the mean squared mass at site k / for Id, sym metric aggregation model for L = 32, 64, 128, 256 
for (a) p = 1.5 and (b) p — 0.5. The numerical fit y = A(p) \J x(l — x) where x = \k\/L is shown in broken line with 
A(p = 1.5) = 1.52 and A(p = 0.5) = 0.34. Set of parameters used: w=l. 



To lowest nonvanishing order in L, the above equations give (Mo) ~ L d and sq < 1 which are the signatures of 
the UA phase. Thus MFT predicts that the system exists only in the UA phase for all p and w (see Fig. ^). Our 
numerical simulations show that this prediction is true at least qualitatively in all cases except in the Id, biased case. 
In Fig. m we show our simulation results for (Mo) in the absence of bias in Id and 2d plotted along with the MFT 
results, and find qualitative agreement. We also measured sq for various system sizes and densities, fixing w = 1 and 
find that it is independent of L in accordance with Eq.(llb). For p = 0.5, we found that SQ UTn ~ 0.42 to be compared 
with s$ IFT ~ 0.33; for p = 1.5, numerically s% um ~ 0.84 whereas s$f FT = 0.6. 



B. UA phase in one dimension 

In this Section, we analyse the UA phase in more detail in the ld, unbiased case. We solve Eq.([3a]) and Eq. ([Sb]) 
for the mean-squared mass in the bulk in this case and find that the MFT seems to violate the conservation law. We 
comment on this limitation of MFT. We also present numerical evidence which shows the simultaneous presence of 
more than one infinite aggregate in the system. 

For convenience, we will choose L to be odd. The mean-squared mass at sites other than the origin obeys the 
following equations, 

2{ml) = (mti) + K+ x ) + 4(m) 2 , \k\ > 2 (12a) 
2{m\) = {ml) + (to) + 4(m) 2 (12b) 
2(to 2 _ 1 ) = (m 2 _ 2 ) + (to) + 4(to) 2 . (12c) 

This set of equations can be solved and one obtains 

(ml) = 2(mf(L-l) + 2(mf(\k\-l)(L-\k\-l) + (m),k^0 . (13) 

To leading order in L, the above result can be written in the scaling form (to 2 .) = L 2 f(\k\/L) where f(x) — x(l — x) 

The mean squared mass (to 2 ,) is actually a measure of the typical mass at site k. To see this, note that we may 
define two types of average over mass distributions at site k, namely, an average (...) over all mass occupations 
including = 0, and an average ((...)) over mass occupations excluding = 0. It is straightforward to see that 
((...)) = (...)/sfc where Sk = 1 — Pk{0) is the probability that site k is occupied. For instance, the average mass (mfc) 
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FIG. 5. Data collapse for the conditional average masses (Mmix) and (Mmix) of the clusters of largest and second largest 
mass in the bulk when the mass at the chipper is Mo for Id, symmetric aggregation model for L — 32, 64, 128, 256. Set of 
parameters used: w=l, p = 1. 



is same on all sites while the typical mass {(rrik)) — (mk) / Sk will depend on k. Similarly, ((m 2 ,)) is the square of 
the typical mass at site fc, implying that (m\) = ((m 2 ,)) may be interpreted as a measure of the typical mass at 
site k. 

We performed numerical simulations to test the mean field prediction for the scaling form (m\) = L 2 f(\k\/L) 
where f(x) = x(l — x). As shown in Fig. ||, we obtain a data collapse with the scaling form (m\) — Lg(\k\/L) with 
g(x) ~ x(\ — x) as the numerical fit for the scaling function. Since g(x) ~ y/x for x close to zero and is constant 
for x 1=3 1/2, the typical mass scales as \L for k close to the origin but as L at sites diametrically opposite to the 
origin. This points to the existence of an aggregate in the bulk as well, consistent with the nomenclature Unpinned 
Aggregate Phase. 

One notes that the MFT solution seems to violate the conservation law since it predicts typical masses of 0(L 2 ) 
at sites situated O(L) away from the chipper. On the other hand, one can also check that Eq.(||) is true using the 
MFT results for (Mo) and sq. Thus MFT seems to be able to describe the vicinity of the chipper more correctly than 
the bulk. The reason for this could be that the mass fluctuations about the mean near the chipper are small due to 
fragmentation unlike those in the bulk where only the aggregation move operates. At a distance of O(L) away from 
the chipper, to a good approximation, one can neglect the presence of the chipper. Then as all the mass resides only 
on one site, the mass-mass correlation function at two different sites (rriimj) is exactly zero in the steady state, in 
strong contrast to the mean field approximation (rriirrij) = (to) 2 . Thus a more refined approximation is required in 
the regions where the aggregation move dominates. 

In Fig. ^, we present the numerical evidence which indicates the simultaneous presence of more than one infinite 
aggregate in the system. We measured the conditional average masses (Mmix) and (Mmix) of the clusters of the 
largest and second largest mass in the bulk respectively, given there is a mass Mq at the chipper. We find that 
data collapse is obtained with the scaling forms (Mmix) = Lf\{Mo/ L) and (Mmix) = Lf2(Ma/L), where the scaling 
functions fi(x) and f2(x) decay almost linearly. Thus a localised infinite aggregate at the chipper and one or more 
mobile infinite aggregates in the bulk can be present at the same time; we would expect more than one infinite 
aggregate to be present in higher dimensions as well. 

In any dimension d, far away from the chipper site, the state resembles that in the absence of the chipper. In the 
latter case, there is a single mobile infinite aggregate which is equally likely to be present at any site with a probability 
1/L d . Although this mobile aggregate with mass of order L d arrives at the chipper infrequently, the probability of 
occupation s of the chipper is of 0(1) - this enhancement occurs because mass can leave only one unit at a time, 
so that it stays for a time of at least order L d . It is implicit in the above argument that an infinite aggregate can be 
formed in the bulk before it hits the chipper. However, this fails to be true in the Id, biased case which explains the 
absence of the UA phase in that case. 
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IV. ASYMMETRIC SINGLE-CHIPPER AGGREGATION MODEL IN ONE DIMENSION 

The mean field prediction that the system exists in the UA phase fails in Id, biased case. In the presence of a 
drift velocity, the one-dimensional system undergoes a phase transition from the NA phase to the PA phase as p is 
increased, keeping w fixed. This exceptional case is the subject of this Section. Although the phase transition survives 
for all nonzero bias, we will only discuss the extreme case when the mass moves only forward (i.e. infinite bias). 

The time required to form an aggregate with mass of 0(L) in the bulk is 0(L 2 ). But in this case, due to ballistic 
motion, the mass clusters return to the chipper in time of O(L) ruling out the formation of an infinite aggregate in 
the bulk. For small w, the sublinear mass arriving at the chipper cannot leave it easily and is temporarily trapped 
giving rise to a localised infinite aggregate (PA phase). As w is increased, the mass leaves more frequently rendering 
the trapping less effective; also this chipped off mass cannot return before the chipper gets empty so that an infinite 
aggregate cannot be sustained at the chipper for large w (NA phase). Thus there is a phase transition in the Id, 
biased case as w (or alternatively p) is varied. 




t 

FIG. 6. Instantaneous mass at the origin Mo(t) vs. t in the steady state of Id, asymmetric aggregation model, (a) The PA 
phase (p = 2). (b) The critical point (p = 0.96875). (c) The NA phase (p = 0.5). Set of parameters used: L—32, w=l. 
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The critical point is determined exactly to be at p c — w in the thermodynamic limit by setting the LHS to be zero 
and so = 1 in Eq.(||). For p < p c , the system exists in the NA phase which is characterised by sq = p/w and (Mo) 
growing sublinearly with L. For p > p c , the system exists in the PA phase in which sq is pinned to its maximum 
value 1 and (Mo) grows linearly with L (see Fig. |J). The probability Pq(0) — 1 — so that the origin is empty serves 
as an order parameter. Our numerical simulations indicate that to lowest order in L, 

( 0{e- L ' La ) , p>w 
I - so = P o (0) ~ I OiL?- 1 ) , p = w (14) 
I constant , p < w 

Since -Po(O) varies continously from the NA phase to PA phase as p increases, this is a second order phase transition. 

Before delving into details, we first present a pictorial representation of the phases. Let us imagine monitoring the 
instantaneous mass at the chipper site in the steady state. Figure ^ shows the time series for the mass at the origin 
as a function of time in the two phases and at the critical point. The mass at the chipper increases due to the mass 
input from the site behind the chipper, and decreases due to fragmentation. Let to and t\ respectively denote the 
consecutive number of time steps during which the chipper is empty and occupied. Then so can be related to to and 
t\ through Sq" 1 — 1 = {to) / {t\) . Depending on whether the ratio (to)/(t\) is zero or not in the thermodynamic limit. 
s is either pinned to 1 (PA phase) or strictly less than 1 (NA phase). In the PA phase, (ti) / (t ) ^> 1 holds due to 
long cascades of successive mass inputs, thus enabling the origin to maintain so equal to 1. The long t\ stretches (see 
Fig. |](a)) enables the chipper to build an aggregate with mass of order of system size. In contrast, in the NA phase, 
(to) ~ (ti) holds as depicted in Fig. ||(b) by the stretches of time during which the origin is empty (due to absence of 
multiple inputs) which reduces so from its maximum value 1. Since t\ is typically not very long, the chipper cannot 
sustain an infinite aggregate in NA phase. Finally, at the critical point, there are multiple mass inputs to the chipper 
but these cascades are not very long as shown in Fig. ^(c). We now turn to a systematic discussion of each phase. 



A. The Pinned Aggregate Phase 

Since the total mass in the system is conserved and the L — 1 sites in the bulk have nonzero average mass (m) and 
sublinear fluctuations about the mean mass (due to absence of a cluster of mass of order L in the bulk as argued above), 
the probability Po(m) that the origin has mass m is ~ \j\f~L exp(— (m — mo) /L) where mo ~ M — (m)L which gives 
Po(O) ~ exp(— L). Thus in this phase, the chipper site is occupied with probability one in the thermodynamic limit 
so that it acts as a reservoir of particles for the bulk. One can think of the system as a semi-infinite, one dimensional 
lattice with a perfect, localised source at the origin injecting monomers into it at the rate w. 

The problem of aggregation in the presence of such a source has been considered previously as well. Some properties 
were studied in |2(J using the technique of interparticle distribution function (IPDF) introduced in and in Q 
by mapping this Id problem to a bounded random walk in 2d. Here we calculate the steady state mass distribution at 
site k denoted by Pk(m) by the generating function method. We define the r-point characteristic function for site k 

at time t, zi- k \x,t) = (e X ^j= k m ^*^). The time evolution equations obeyed by Z^{X,t) with a perfect, localised 
source at k = are given by 



' Z ' = Z%\ 1) + Z {k \ - 2Z« , k ? 1 , r + (15a) 



Of 



dZ'r 



(1) 



Ot 



= Z^\ + {we- x - 1 - w)Z^ , r^O (15b) 



with the boundary condition Zq (A,i) = 1 for all k > 0. In the steady state, zi k \\,t) is independent of t. 

We need to solve for z[ (X) which is the Laplace transform of Pfc(m) with respect to m. One can easily solve 



for Z^iX) using Eq.(|15bD and the boundary condition Zq = 1 which on taking the inverse Laplace transform for 



(i) 



large m gives P\(m) = e m ' w /w. Thus the mass distribution at the first site decays exponentially. Now to find the 
probability distribution in the bulk, we define 

H(x,y,\)= E Zi k HX)x k y r . (16) 

k=2,... r=l,... 
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FIG. 7. Data collapse for the unnormalised probability distribution of mass at site k = 64 (solid squares), 128 (crosses), 256 
(squares) in the PA phase plotted against the analytical result. Inset: log-log plot of probability distribution for nonzero mass 
at site k plotted against the analytical result. Set of parameters used: L = 256, p = 2, w = 1. 



Using equations Eq. (|15a| ) and Eq.([15 

H[x,y, A) 



I, one obtains 

xy 



y — 2y 



G(x, A) — xy 



1-x 



i-fy 



(17) 



where 1// = 1 + w — we " A and G(x, X) = £ fc=2 Z i 1 k) (X)x k . Clearly, the inverse Laplace transform of G(x, X) w.r.t. 



A and x gives Pk(m) for k > 1. To calculate G(x, A), we use the same method as in 23 1. We first note that the 
denominator of H (x, y, X) has roots at y± = 1 ± \/l — x. The root at y — y- lies inside the unit circle which in 
the real space (i.e. r-space) gives an exponentially increasing solution which is disallowed since the inverse Laplace 
transform of H(x, y, A) gives a probability, which is always bounded above. To avoid this pathological solution, we 
demand that y = j/_ is a zero of the numerator as well, i.e. 



G{x, X) ~ xy- 



1 



(18) 



1-x ■ l-fy- 

For A — > and x — * 1, the inverse Laplace transform can be easily found. We find that Pk(m) has a scaling form 
4>{u)/m 2 where u = m/2wy/k and the scaling function is 

Aw 



4>{u) 



-.u 3 e- u 



(19) 



One can also obtain the probability distribution Pfc(O) of zero mass at site k by solving for G(x, X) when A — ► oo 
and x — > 1. To leading order in k, one obtains Pk(0) = 1 — Sk = 1 — l/Vftk. We tested our calculations against 
the numerical simulations performed on the single-chipper model and found reasonably good agreement (see Fig. |^). 
Thus the typical mass at site k grows as 

We now calculate the typical spacing between the masses using the method of IPDF (22| . We define Ek {r, t) as the 
probability that the sites k to k + r are empty (including k and k + r). Then Ek(r, t) satisfies the following equations, 



dE k (r,t) 
dt 

dE 1 {r,t) 
dt 



= E k (r-l,t)-2E k (r,t)+E k _ 1 (r + l,t) , fc^l 
= E 1 {r-l,t)-(l + w)E 1 (r,t) , 



(20a) 
(20b) 



with the boundary condition Ek{— = 1 for all k > 0. In the steady state, Ek(r,t) is independent of t and one 
can write the above equation in terms of F±(z,y) = J2k=2 Sr=o Ek(r)z k y r and R(z,r) = J2k=2 E k (r)z k , 
r = 0,1 We find 
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Fi(z,y) 



y 

y 2 - 2y + z 



zR(z,0) 



(21) 



y (1 + w) (1 + w — y) 1 — z 

With the same reasoning as in the previous calculation, we demand that the numerator of Fi(z,y) evaluated at 
y_ = 1 — yl — z be zero. This condition gives 

1 1 



R(z,0) 



l-z (1 + w)(l + w - y-) 



z y- 



(22) 



which on inverting the Laplace transform gives the probability Ek(0) which is same as P k (0). One can check that the 
result for -Pfc(O) quoted above is reproduced. Substituting R(z, 0) in Fi(z, y), we obtain 



Fi(z,y) 



y - y+ 



1 — z (1 + w — y)(l + w — y-) 



(23) 



where y + = 1 + Vl + z. 

We further define D k (r) as the probability that both sites k and k + r are occupied but no sites in between. Then 
D k (r)=E k+1 (r-2)-E k+1 (r-l)-E k (r-l) + E k (r) , fc^l.r^O . (24) 
Defining the Laplace transform of D k (r) with respect to k and r as F 2 (z, y) = ^2 k=2 X)r=o D k (r)z k y r , we obtain 



F 2 (z,y) = ~[R{z,-2 



+z(l-y) 1 + 



yz 

1- z 

y 

2- 2/ 



z(l + z) 
1 - z 



+ F 1 {z,y)[ y - + (l-y)(l + ±) 



(25) 



(1 + tu) (1+^-2/) 



where F\{z,y) was calculated above. 

Then the typical empty space (r k ) in front of site k can be obtained by taking the inverse Laplace transform of 
(dF2/dy)\ y= i with respect to z. One finds 



we w2k eric(wVk) 



(26) 



which for large k implies that the typical empty space ((rk)) — (r k )/s k in front of occupied site k varies as 

Thus, in the PA phase, the mass in the bulk is distributed in \fL clusters with typical mass \fk at site k, and an 
empty stretch of length \fk in front of it. 



B. The Non Aggregate Phase 

As discussed earlier, in the NA phase, sq is strictly less than 1 due to substantial time stretches of typical length 
(to) during which the origin is not occupied (Fig. ||). Thus the origin does not act as a perfect source unlike in the 
PA phase and we could not calculate the mass distribution in this phase. However, one can obtain useful information 
about the nature of this phase by simple scaling arguments. We begin by arguing that (Mo) grows sublinearly with 
L and then provide numerical evidence for the absence of an infinite aggregate (i.e. mass proportional to L) in the 
bulk, thus concluding that there is no infinite aggregate anywhere. 

Monte Carlo simulations indicate that (M ), (to) and (ti) vary as a power law in L with a constant as a next order 
correction (see Fig. ||). Therefore, we assume (M ) = aL 13 + b, (t ) = cqL 1 + d and (tx) = ciL 11 + d\. Solving for s 
using Eq.(||), we obtain 

•o = ~. + O (t^s) ■ (27) 



w 



Then since the lowest order term for so is a constant and s 1 — 1 = (to)/(t\), one obtains 7 = rj. Using this identity 
and retaining terms to lowest order in L we further obtain 

*o = -^— + o(±-) . (28) 



c + ci V^ 7 



Comparison of Eq.(27) and Eq. (B8h yields + 7 = 1. 
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FIG. 8. Data collapse for the typical mass at site k 7^ in the NA phase for L — 64, 128, 256, 512, 1024. Set of parameters 
used: p = 1, w = 2. In this figure, all the sites are labelled by positive integers in the clockwise direction with the chipper at 
the origin. 

Now we consider the typical mass ((nik)} in the bulk where the site index k is labelled by positive integers in 
the clockwise direction with the chipper at the origin. Numerically, ((mfc)) is observed to obey the scaling form, 
((?7ifc)) = Vkf (k/L) as shown in the data collapse in Fig. |§|. The scaling function is a slowly varying function, which 
gives sublinear mass everywhere in the bulk. Thus in the NA phase, the mass is distributed in clusters of typical mass 
~ 0(VL) for k ~ 0(L). 

Since the mass at the site just behind the chipper is of 0(\/L) and as one can see in Fig. [| the number of mass inputs 
to the chipper is typically of O(l), it follows that (ti) ~ 0(V~L) so that r\ = 0.5, which further gives /3 = 7 = r\ = 0.5 
due to the two scaling relations above. These exponent values are checked numerically, as shown in Fig. H. 
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FIG. 9. Plot of (Mo), (to) and (ti) as a function of \[L in the NA phase. Since all the three curves are linear, it follows that 
(Mo) = aL 13 + b, (to) — cqL 1 + do and (ti) — c±L v + d\ with /3 — 7 = r\ — 0.5. Set of parameters used: p = 1, w = 2. 
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C. The Critical Point 



We studied the critical point mainly numerically by studying the L dependence of (M ), (t ) and (t\). Assuming 
a power law dependence for (Mo) ~ L 13 , (to) ~ L 1 and (ti) ~ L v , a naive best fit in the range L = 32 to 2048 gives 
(3 ~ 0.62, 7 ~ 0.38 and 77 ~ 0.80. However, the effective exponent calculated using the successive ratios of L used in 
the simulations shows systematic trends: (3 decreases as L increases while both 7 and r\ increase with L. 



0.56 



0.54 



0.52 



0.5 




1000 



1500 



2000 



2500 



FIG. 10. Plot of (M )/L - 62 (squares) and (M )/VZ (log L) 13 ' , ft ~ 0.7 (triangles) as a function of L at the critical point 
in the lei, biased case. The variation of (Mo) as \fL with multiplicative logarithmic correction seems to be a better fit than 
the fit to a pure power law. The plot of (Mo)/L ' 62 vs. L has been scaled by a constant factor. Set of parameters used: 
p=l,w = pL/(L-l). 
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FIG. 11. Plot of (to)/L°' i8 (squares) and (to) (log L) 1 /V~L, 7' ~ 0.74 (triangles) as a function of L at the critical point 
in the lei, biased case. The variation of {to} as yl with multiplicative logarithmic correction seems to be a better fit than 
the fit to a pure power law. The plot of {to)/L°' S8 vs. L has been scaled by a constant factor. Set of parameters used: 
p=l,w = pL/(L-l). 
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A good fit is obtained if one allows for logarithmic corrections in the power laws for (Mo) and (to) with (3 = ^ = 1/2 (see 
Fig. m and Fig. |l|). This indicates that both (M ) and (t ) vary as \[L with multiplicative logarithmic corrections. 
Using Eq.(||) and the condition of criticality, one can easily show that 77 + j3 — 7 = 1. Using /3 = 7 = l/2as suggested 
by the discussion above, we obtain 77 = 1. The above quoted values for the exponents f3 and 7 are consistent with the 
scaling relation (3 + 7 = 1 which seems to be true in the NA and the PA phase as well. In the NA phase, we have seen 
in Section IV B that this scaling relation holds. In the PA phase, we found numerically that (to) is independent of L 
in the thermodynamic limit implying that 7 = 0; further since (3 = 1, it follows that the scaling relation (3 + 7 = 1 
holds in this phase also. 

We have also numerically studied P{t\) which is the probability that the origin is occupied for t\ consecutive time 
steps. At the critical point, this distribution shows two peaks - a broad peak occurs at the typical time scale t\ ~ L 13 
while there is a narrow peak at t\ — N p /w. We have not been able to reliably separate out these two contributions 
to P(t\) but since the scaling relation yields rj = 1, one may expect that the second peak dominates the asymptotic 
value of 77. 



V. SINGLE-CHIPPER FREE PARTICLE MODEL 



A. The Model 



In this Section, we study the steady state of a model of noninteracting particles in which a localised infinite aggregate 
is formed solely due to disorder. In this model the particles diffuse freely in the bulk except at certain, quenched sites 
referred to as the chipper sites. It can be shown that this system exhibits a phase transition from the NA phase to 



the PA phase (as defined in Section II B ) as the particle density is increased in all dimensions for all bias. In this 
Section, we will demonstrate this result in the presence of a single chipper. We will see that even a single defect is 
capable of inducing a phase transition in all dimensions, unlike in the aggregation model. We will discuss the solution 
with an extensive number of chipper sites in Section [v|. 

Our model is defined on a e?-dimensional hypercubic lattice with periodic boundary conditions on which we consider 
the biased and unbiased diffusion of a conserved number of particles. At any site except the chipper, each particle 
attempts to hop out at rate 1. Since the particles do not interact, the hopping rate out of site k occupied by m 
particles is m. At the chipper site, assumed to be at the origin, the hopping rate is a constant w, independent of the 
number of particles. Thus, the rate Uk(m) at which a particle leaves site k which has m particles is given by 

Uk(m) = w 4,o + m(l - 8 k, a) , m^O . (29) 

In the fully asymmetric, Id case, the model described above can be mapped onto a traffic model on a one lane road 
with no overtaking allowed. We represent each site in this model as a car and each particle as a vacant site. Then 
a system of M free particles on a lattice of size L maps onto a system of L particles with hard core interactions on 
a lattice of size L + M. In this new representation, the special car (corresponding to the chipper site) moves with a 
constant rate, irrespective of the headway in front of it and rest of the cars (sites other than the chipper) move with a 
rate proportional to the headway in front of it. Note that the sitewise disorder in the free particle model corresponds 
to particlewise disorder in the traffic model. 



B. Phase Transition in arbitrary dimensions 

The steady state of this system can be found exactly by noting that the hopping rates Eq.(^9|) in this model 
correspond to a special choice of rates in the Zero Range Process (l^Jl^]. In this process, a particle at site k hops 
to its nearest neighbour independent of the state at the target site, so that the interaction has zero range. The 
steady state measure of this process can be found in any dimension with or without bias. A convenient way to find 
it is by using the condition of pairwise balance J24|] which states that for any given configuration C, one can find a 
configuration C in one-to-one correspondence with C" such that 

W(C -> C) P(C) = W(C -> C") P{C) , (30) 

where P(C) is the probability of a configuration C in a system with N sites and W(C — » C) is the transition rate 
from C to C The above condition is satisfied by 
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1 N 

P(C) = Tr l[fk(m) , (31) 



«k- 



where 



m 1 

/«=n^) for m >° 

= 1 for m = , (32) 

and Ufe(m) is the rate at which a particle hops out of a site fc having m particles and Af is the normalisation constant. 
fk(jn) is defined upto a multiplicative factor v m where v can be interpreted as the fugacity. 



For our model, using the above solution for the particle distribution /fc(m) and Eq.(29), we find that the normalised 



probability distribution Pfc(m) that the site k in the bulk has m particles has the Poisson form, 

e~ v v m 

P k (m) = P{m) = ^,fc^0 . (33) 

m! 

The fugacity v will be determined using the conservation law. At the chipper site, one can obtain solution for Po(m) 
when v < w, 

prM = (i--)(-) m , -<i , (34) 

where NA in the superscript stands for the Non- Aggregate phase which we have added in the anticipation of a phase 
transition. 



Using equations Eq. (|33[) and Eq.(34), one can easily see that the average particle number at site k ^ is (mk) 
v = wsq which further leads to constraint equa tion Eq.Q. Since this constraint equation is identical to that in the 



aggregation model, then as discussed in Section IIE , the steady state of this model also has three possible phases (see 
Fig. ||). But, as we will see, this system never exists in the UA phase and there is a phase transition from the NA 
phase to the PA phase as p is increased, keeping w fixed, in all dimensions, irrespective of bias. 

Solving for the average number of particles at the origin (Mo) using Eq.(34) and substituting in Eq.([s]), we obtain 



so 



1 - s 



pL d - ws (L d - 1) , < s < 1 • (35) 



To leading order in L, one obtains two solutions, namely, sq — p/w and sq = 1. Since the above analysis holds only 
when so is strictly less than 1, the only valid solution is so = p/w for p < w. For p > w, sq is pinned to its maximum 
value 1. Thus there is a phase transition from the NA phase to the PA phase at p c = w, as p is increased, keeping w 
fixed. 

One can solve for Po(m) in the PA phase and at the critical point using the conservation law Eq.(||) and the 
fact that Eq.(|33|) is valid for all so in the range [0,1]. Since the total number of particles is conserved, Po{m) = 
^2 C P(C) S(J2k^o mk = M — m). The quantity on the RHS can be calculated strightforwardly. The result is 



pa, s / 1 ( (M -m- (m)(L d -l)f , 

p ° ^ = v ^ a )(^-D exp ( ( ^-i)(U 1 ■ 



where (m)=v and (m 2 )~v(l + v) from Eq. (p3|) . One can see from the above equation that (Mo) = M — w(L d — 1) ~ L d 
and Pq >a (0) = 1 — so ~ e~ L in the PA phase. Further, at the critical point, we have pL d = w(L d — 1) using which 
in Eq.(|36]), one obtains Po(m) at the critical point, 



P " {m) ^ V n(mi)(L d -l) CXP [-{L*-l)(m*) ) ' (37) 

which gives a power law decay in L for Pq P (0). 

The phase transition is brought about in this model in the following way. In the NA phase which is a low density 
phase, the typical number of particles at all sites including the chipper is of 0(1). As the density is increased, there 
is a phase transition to the PA phase at p c — w. In this high density phase, each site in the bulk still supports only 
0(1) number of particles in accordance with Eq.(|33"|) so that the extra particles condense on the chipper giving rise 
to (Mo) ~ L d . The mechanism of phase transition in this model is similar to that in Bose-Einstein condensation as 
was pointed out in M] in a similar lei model which shows a phase transition with Ufe(m) = it? 5k,o + (1 — <5fc,o)- 
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VI. EXTENSIVE DISORDER 



In this Section, we describe a possible scenario in the more interesting and physically relevant situation when there is 
an extensive number of chipper sites. These are assumed to be placed randomly, with quenched random chipping rates 
Wk at site k distributed according to a distribution Prob(iUfc). We are interested primarily in the aggregation model 
but as argued below, on large space and time scales the behaviour of this model with extensive disorder resembles 
that of the corresponding free particle model (a generalization of the model in Section ^ with an extensive number 
of chippers) which is exactly solvable. 

Let us consider the diffusion-aggregation process with no bias and take the initial condition to have a random 
distribution of masses. Then the finite concentration c of chipper sites brings in new length and time scales into 
the problem, namely the mean spacing l c ~ c 1 ^ between the chippers and the associated diffusion time t c = 1%. 
Evidently, l c and t c define respectively the relevant length and time scales over which the diffusing mass clusters sense 
the presence of an extensive number of defects. 

Let us consider a finite but low concentration of chipper sites so that 1 < l c < I. On time scales t <C f C) we would 
expect the system to behave roughly as independent, finite systems of size l c with typically a single c hipp er each. The 



typical state on these short time scales thus resembles the UA-like steady state discussed in Section |IlJj. The typical 
mass of clusters, both mobile and localised is limited by the time and grows proportional to t d l 2 as long as t -C t c . 

As the time crosses t c , the finite spacing between the chippers becomes relevant. The primary effect is to limit the 
size of the aggregates formed in the non-chipper region by the diffusion-aggregation process to ~ if, since on a time 
scale > t c the mass cluster is likely to encounter a chipper and get trapped. For t ^> t c , the coarse grained view of 
the system is that of mass exchanged between close- by random rate chippers with each exchange taking a time ~ t c . 
To the extent that only finite aggregates (with mass ~ if <C L d ) are formed in the transit between chippers, it is 
plausible that on large length and time scales, we may ignore interaction effects (i.e. coalescence) and think of the 
system as effectively free-particle like in the non-chipper region. 

The free particle model with extensive disorder is solvable along the same lines as the single-chipper problem 
described in Section |y|. It defines a Zero Range Process with the hopping rates, 

Wfe(m) = Wk , if k is a chipper site (38) 
= m , if k is not a chipper site 

Let x denote the fraction of chipper sites. We recover the single-chipper model as x — > 0, while x — > 1 corresponds 
to every site being a chipper site, and is the model considered in ]9|,|l0|]. For all x, in the low-density phase, the mass 
distribution on the non-chipper sites follows the probability distribution of Eq.(|3^), while at chipper sites, Eq.(|3~4|) is 
valid with w replaced by Wk- Let Sk be the occupation probability of site k, and so refer to the occupation probability 
of the site with the lowest chipping rate, wq — Min{wk}- Then in the steady state, the spatial uniformity of the 
current leads to (i) w^Sk — wqSq if k is a chipper site, and (ii) (mk) = woSq if k is not a chipper site. Using these 
relations, one can write the mass conservation equation analogous to Eq.(|35|) as 

1 s , f , Prob(w)s 

' aw — h (1 — x)s W() — p . (39) 



L d 1 — so J w/wq — so 

where we have separated out the first term corresponding to the slowest chipper. 

One can analyse Eq. (|39"| ) for various x as follows: (a) In the limit of a single chipper (reached as x approaches 0), 
we know from Section M that as the density is increased, there is a phase transition from the NA phase to the PA 
phase with an infinite aggregate at the chipper site, (b) In the limit x — > 1, the model reduces to that considered in 
P,[To|, where it is shown that for Prob(u;fc) ~ (wk — Wo) u as Wk — > wq : the system stays in the NA phase for all p if 
v < 0. On the other hand, if v > 0, then there is a transition to the PA phase with an infinite aggregate at the site 
with chipping rate wo when the density crosses the critical density given by 

Pc (x = l)= [dw^l . (40) 

J W/Wq-I 

(c) For < x < 1, there is a transition from the NA to the PA phase if Probfui) is chosen as in (b) above. The critical 
density can be determined by taking so — > 1 and L d (l — sq) — > oo in Eq . (P9[) , with the result 

p c (x) = xp c (x = 1) + (1 - x)w . (41) 

Thus the critical density interpolates linearly between its values in the limits x = and x = 1. 
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In view of the correspondence discussed at the beginning of this Section, we would expect the aggregation model 
with a similar distribution of chipper sites to show a phase transition from the NA to the PA phase. In the former 
phase, there are both localized and mobile aggregates with typical mass ~ if. In the PA phase, the distribution of 
masses is similar to that in the NA phase at all sites except the slowest chipper with Wk = wq\ at this slowest site, 
there is an aggregate with mass of order volume. 



VII. SUMMARY 



In this paper, we introduced a minimal model to study the effect of quenched, sitewise disorder in a aggregation- 
fragmentation system. Our model had some simplifications: the fragmentation was allowed to occur only at the 
trapping sites, and mass-independent kernels for aggregation and fragmentation were considered. Despite these 
simplifications, it retains the important physical effects of diffusion, aggregation, fragmentation and trapping. 

We studied the case of a single-chipper aggregation model in detail. In all cases except the Id, biased case, the 
system exists in the UA phase in which an infinite aggregate localised at the chipper site can coexist with mobile 
infinite aggregates in the bulk. The simultaneous existence of more than one infinite aggregate is a new feature absent 
in previous studies of translationally invariant systems. In the Id, biased case, there is a phase transition from a phase 
in which a localised aggregate is formed at the chipper site (PA phase) to the one in which no aggregate is formed 
anywhere in the system (NA phase) as the density is decreased. 

We also studied a variant of the above aggregation model in which particles chip off at a single site but diffuse freely 
in the bulk. This model can be solved exactly and shows a phase transition from the PA phase to the NA phase as 
density is increased in all dimensions for all bias. 

Finally we discussed a likely scenario for the aggregation model in the presence of extensive disorder and argued 
that interaction effects arising due to coalescence can be ignored on large enough time scales. It would be interesting 
to check this expectation by a more detailed study of the system with extensive disorder. 
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